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^ ■ The ISM of active galaxy centers is exposed to a combination of cosmic ray, 

FUV and X-ray radiation. We apply PDR models to this ISM with both 'nor- 
mal' and highly elevated (5 x 10~ 15 s _1 ) cosmic-ray rates and compare the results 
^ ■ to those obtained for XDRs. Our existing PDR-XDR code is used to construct 

models over a 10 3 — 10 5 cm -3 density range and for 0.16-160 erg s _1 cm -2 im- 
pingent fluxes. We obtain larger high J (J > 10) CO ratios in PDRs when we 
use the highly elevated cosmic ray rate, but these are always exceeded by the 
corresponding XDR ratios. The [CI] 609 /zm/ 13 CO(2-l) line ratio is boosted by 
a factor of a few in PDRs with n ~ 10 3 cm -3 exposed to a high cosmic ray 
rate. At higher densities ratios become identical irrespective of cosmic ray flux, 
while XDRs always show elevated [CI] emission per CO column. The HCN/CO 
and HCN/HCO + line ratios, combined with high J CO emission lines, are good 
diagnostics to distinguish between PDRs under either low or high cosmic ray ir- 
radiation conditions, and XDRs. Hence, the HIFI instrument on Herschel, which 
can detect these CO lines, will be crucial in the study of active galaxies. 



Subject headings: ISM: X-rays — ISM: Cosmic Rays 
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1. Introduction 

In centers of late-type galaxies, such as M 82, NGC 253, and Maffei 2, molecular line 
intensity ratios 1 are frequently found to be high, e.g. CO(4-3)/CO(l-0)~ 60 (Israel et al. 
1995; White et al. 1994; Israel & Baas 2003) requiring high gas densities n > 10 5 5 cm -3 , and 
high temperatures of T > 50 K. However, FUV photons are easily attenuated by dust and 
do not penetrate very deep into clouds and the galaxies are observed with beams covering 
regions typically hundreds of parsecs in size. At such large spatial scales, excited dense gas is 
not expected to have very large filling factors. UV radiation seems incapable of maintaining 
very large gas fractions at temperatures of T ~ 50 — 150 K, and regular PDR models do not 
explain the observed high ratios. In addition, Israel & Baas (2002) have observed the [CI] 
609 fim line in the centers of late-type galaxies, and measured [CI] 609 /zm/ 13 CO(2-l) line 
intensity ratios in the range 20 - 60, which are hard to explain by (low-CR) PDR models. 
The C + /C/CO transition zones in these PDRs contain only a thin layer in which neutral 
carbon has a high abundance. At the cloud edge, carbon is ionized, and deep in the cloud 
all carbon is locked up in CO. 

Various authors have invoked elevated cosmic ray fluxes caused by greatly enhanced 
supernova rates in galaxy centers in order to explain large molecular gas masses at high 
temperatures (Suchkov et al. 1993; Bradford et al. 2003), and to explain high [CI] intensities 
and column densities relative to 12 CO and 13 CO (Pineau des Forets et al. 1992; Schilke et al. 
1993; Flower et al. 1994). As most cosmic rays are produced in supernovae, their flux is 
proportional to the star formation rate, which is about 1 M Q yr" 1 for the Milky Way. In 
circumnuclear starbursts, star formation rates may be two orders of magnitude higher or 
more. Such galaxy centers may also contain an embedded accreting black hole producing 
X-rays. Like cosmic rays, but unlike UV photons, these X-rays can also penetrate through 
large column densities (N H > 10 24 cm -2 ), and can cause the observed high line ratios over 
areas as large as 500 pc, when the emitted flux is high enough (Meijerink & Spaans 2005; 
Meijerink et al. 2006). 

In this paper, we investigate whether in galaxy central regions PDRs with and without 
enhanced cosmic ray fluxes can be distinguished from XDRs on the basis of observable 
atomic and molecular line ratios. To this end, we calculate line intensities for PDRs with 
very different cosmic ray rates and compare the results to those obtained for XDRs with the 
same radiation fields and column densities. 



1 Intensity ratios of lines a to lines b are related to brightness temperature ratios by the cube of the line 
frequencies: {v a /vb) 3 
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2. PDR and XDR models 

We have constructed a set of PDR and XDR models from the codes described by Mei- 
jerink & Spaans (2005) and Meijerink et al. (2006), in which we varied both the incident 
radiation field and the density. The thermal balance (with line transfer) is calculated self- 
consistently with the chemical balance through iteration. Absorption cross sections for X- 
rays (1-100 keV) are smaller, ~ 1/E 3 , than for FUV photons. Therefore, PDRs show a 
stratified structure while the changes in the chemical and thermal structure in XDRs are 
very gradual. Species like C + , C and CO co-exist in XDRs, and large columns of neutral 
carbon (unlike in PDRs) are produced. In the XDRs, additional reactions for fast elec- 
trons that ionize, excite and heat the gas are included. The heating efficiency in XDRs is 
much higher. Since we focus on galaxy centers, we have assumed the metallicity to be twice 
Solar. We take the abundance of carbon to be equal to that of oxygen, since the carbon 
abundance increases faster than oxygen for larger metallicity. The precise C:0 ratio does 
not affect our general results. We have calculated PDR models for both a 'normal' (low: 
( = 5x 10" 17 s" 1 - cf. van der Tak k van Dishoeck 2000) and a high (( = 5 x 10" 15 s" 1 ) 
cosmic ray flux corresponding to a star formation rate of ~ 100 M yr" 1 . The model input 
parameters are summarized in Table 1. The range of free parameters (n = 10 3 — 10 5 cm" 3 , 
Go = 10 2 — 10 5 /Fx = 0.16 — 160 erg s" 1 cm" 2 ) is representative for the conditions in galaxy 
centers. That is, from low-J CO (critical density 10 3 — 10 4 cm" 3 ) and HCN 1-0 (critical 
density ~ 10 5 cm" 3 ) observations it is apparent that gas in galaxy centers must exhibit the 
density range that we model. The irradiation conditions are typical for Milky Way PDRs 
like the Orion Bar (a high mass star-forming region) as well as a generic 10 44 erg s~l Seyfert 
nucleus X-ray luminosity for distances of about 100 pc and up. 

3. Chemical and thermal structure 

Higher cosmic ray (CR) ionization rates do not much affect the chemistry at the cloud 
edge, but large effects occur beyond the H/H 2 transition (see e.g. Fig. 1). As the CR flux 
in the PDRs is enhanced, electron, carbon and hydrogen abundances decrease much less 
beyond the H/H 2 transition, due to larger ionization/dissociation rates. Ion abundances 
also remain higher, causing H 2 and OH to have higher abundances as well. As cosmic ray 
ionization contributes to the gas heating, higher incident CR fluxes raise gas temperatures 
deep in the cloud, with roughly T ~ ( 1 ^ 3 at low densities (n ~ 10 3 cm" 3 ). At high densities, 
n = 10 5 cm" 3 , dust acts as an effective coolant if T > T^, and the rise in the kinetic 
temperature is less pronounced. Both temperatures and abundances in a PDR become 
higher when CR rates are increased and ratios of emergent line emissions are also modified. 
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4. CO line intensities and ratios 

a. Density n = 10 3 cm -3 (Table 2). In high-CR PDR clouds, a large fraction of all CO is 
dissociated, and CO is a factor ~ 100 less abundant than in low-CR PDRs. However, in the 
latter, the CO lines are optically thick and CO line intensities are comparible for the same 
incident flux. In high-CR PDRs, the low-J CO line ratios are somewhat larger than those 
in low-CR PDRs. The higher transitions show more of a difference and are diagnostically 
more valuable. In the XDRs, the CO gas temperature is on average much higher than 
in the PDRs, but marginally higher CO intensities occur only for the weak radiation field 
F x = 1.6 ergs s _1 cm" 2 . In XDRs with stronger radiation fields, CO line intensities are 
much lower than in PDRs, because column densities are too small to attenuate the X-rays 
significantly, and CO abundances are very low throughout the cloud. The very high CO (2- 
1)/CO(1-0) and CO(4-3)/CO(l-0) ratios in XDRs, which can be much higher than in the 
corresponding PDR cases, merely reflect the weakness of the lower J lines. 

b. Density n = 10 4 cm" 3 (Table 3). High-CR PDRs have more intense CO lines than 
low-CR PDRs, but their high- J /low- J CO line ratios are only marginally higher than those 
of the corresponding low-CR PDRs. Again we find that XDR CO line intensities exceed 
those in PDRs in relatively weak radiation fields (G = 10 3 or F x = 1.6 ergs s" 1 cm" 2 ). In 
XDRs with stronger radiation fields, the lower J CO lines (up to CO (4-3)) are much weaker 
than in PDRs, but the higher rotational lines are stronger. All XDR CO ratios are (much) 
larger than even the high-CR PDR CO ratios, and ratios at the intermediate and high J 
levels above CO (6-5) are diagnostically particularly meaningful. 

c. Density n = 10 5 cm" 3 (Table 4). There are no longer significant differences between 
low-CR and high-CR PDRs. However, the XDR CO lines are strong in the lower transitions, 
and even more so at the intermediate and higher rotational transitions above CO (4-3). Their 
ratios are always larger than the corresponding ratios in any PDR model (see Meijerink et al. 
2006 for a more detailed discussion). Some high- J CO intensities and ratios are left blank, 
since no significant emission was found for these lines due to the low fractional abundances 
of CO and the high critical densities of the transitions. 

5. [CI] 609 /im/ 13 CO(2-l) ratios 

XDRs, quite unlike PDRs, have significant neutral carbon abundances; their [CI] inten- 
sities behave as volume tracers. At the same time, XDRs generally have weaker low J CO 
and 13 CO lines than PDRs, for n < 10 4 cm" 3 and Fx > 1 ergs s _1 cm" 2 . Thus, in XDRs 
[CI] 609 /xm/ 13 CO(2-l) line intensity ratios are much larger than in corresponding PDRs. 
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Only at low densities, n < 10 3 cm 3 , do high-CR PDRs behave in a fashion intermediate 
between XDRs and low-CR PDRs. In such low-density, high-CR PDRs, CO and 13 CO dis- 
sociation causes simultaneous 13 CO line weakening and [CI] line strengthening, resulting in 
[CI]/ 13 CO line intensity ratios four times higher than seen in low CR-PDRs, but still (much) 
lower than those seen in XDRs. Even at higher densities (n > 10 4 cm -3 ) and stronger radi- 
ation fields (Go = 10 5 ), the high-CR PDRs fail to produce [CI]/ 13 CO ratios more than 1.5 
times those of low-CR PDRs. However, in XDRs the [CI]/ 13 CO ratios remain very much 
larger and thus provide an excellent tool to distinguish between PDRs and XDRs. 



6. HCN/CO and HCN/HCO+ ratios 

Meijerink et al. (2006) found that HCN/CO and HCN/HCO+ line intensity ratios dis- 
tinguish between 'normal' (low-CR) PDRs and XDRs. These ratios are slightly different in 
high-CR PDRs. At densities of n = 10 4 cm~ 3 , the HCN/CO ratios are lower in high-CR 
PDRs than in low-CR PDRs, especially for the J = 4 — 3 transition, but they become more or 
less identical at densities of n — 10 5 cm -3 . In XDRs, the HCN/CO ratios are almost invari- 
ably significantly lower than in either PDR. High-CR PDRs with densities of n — 10 4 cm -3 
have HCN/HCO + ratios higher than low-CR PDRs in the J = 1 — transition and lower 
than low-CR PDRs in the J = 4 — 3 transition. This opposite behaviour reflects a cosmic 
ray induced shift in the HCN and HCO + abundances to larger (cooler) columns. However, 
at higher densities (n = 10 5 cm -3 ), high-CR PDRs have HCN/HCO + ratios always lower 
than those in low-CR PDRs, mainly because of a boost in the HCO + production. In all cases 
where XDRs produce HCN and HCO + emission observable at all, the HCN/HCO + ratios 
are (much) lower in the XDR than in either PDR. 



7. Conclusions 

1. CO line intensity ratios increase when cosmic ray ionization rates are enhanced, but 
they remain smaller than those in XDRs. In particular high- J (J > 10) CO lines 
(which will become observable with HIFI in ESA's Herschel space observatory) allow 
to distinguish between (high-CR) PDRs and XDRs. Using the HIFI time estimator 
and a beam filling factor of 0.05, we find that a line intensity of 9 x 10 -6 erg s" 1 cm -2 , 
the largest CO(16-15) intensity produced by our PDR models, will be detectable with 
HIFI in about 4 hours, while the very bright lines received from highly irradiated XDRs 
are detectable within minutes. 



2. [CI] 609 /iin/ 13 CO(2-l) ratios are much higher in high-CR PDRs than in low-CR PDRs 
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at modest densities of n — 10 3 cm" 3 . At higher densities of n > 10 4 cm -3 , this 
difference vanishes. In XDRs, the ratios are always larger than in PDRs at the same 
density, independent of incident radiation field. 

3. HCN/CO and HCN/HCO + line ratios are good diagnostics to distinguish between 
PDRs and XDRs. As the ratios obtained for low-CR and high-CR PDRs are different, 
combination with high- J CO lines is both crucial and profitable in the study of the 
active galaxy centers. 

4. Although our model results broadly distinguish between low and high CR PDRs and 
XDRs, there is some degeneracy when constraining the four parameters (density, CR 
rate, Go and Fx) through individual ratios. Therefore, a combination of various ratios 
should be considered. 



We thank J. le Bourlot and M. Elitzur for useful discussions on X-ray versus cosmic-ray 
irradiation, M. Hogerheijde and P. van der Werf for discussions on HIFI sensitivities, and an 
anonymous referee for his constructive comments. 
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Table 1: PDR and XDR models 



Density 


PDR 1 (Go) 


XDR (F x ) 


Size (pc) 


10 3 


10 2 , 10 3 , 10 4 


0.16, 1.6, 16 


10 


10 4 


10 3 , 10 4 , 10 5 


1.6, 16, 160 


1 


10 5 


10 3 , 10 4 , 10 5 


1.6, 16, 160 


1 



1 Both low and high CR PDRs are calculated. 
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Table 2: Line intensities (erg s 1 cm 2 sr x ) and ratios at density n = 10 3 cm 



Model 




Low-CR PDR 






High-CR PDR 






XDR 




Radiation field 


G = 10 2 


G = 10 3 


Go = 10 4 


G = 10 2 


G = 10 3 


Go = 10 4 


F x = 0.16 


F x = 1.6 


F x = 16 


CO(l-O) 


8.2(-8) 


8.6(-8) 


1.0(-7) 


8.4(-8) 


8.3(-8) 


8.1(-8) 


1.2(-7) 


3.7(-10) 


3.0(-10) 


CO(2-l) 


4.8(-7) 


5.1(-7) 


6.6(-7) 


6.3(-7) 


6.3(-7) 


6.4(-7) 


1.0(-6) 


9.3(-9) 


9.1(-9) 


CO(3-2) 


9.4(-7) 


9.9(-7) 


1.5(-6) 


1.6(-6) 


1.6(-6) 


1.6(-6) 


2.8(-6) 


3.4(-8) 


5.2(-8) 


CO(4-3) 


1.4(-6) 


1.6(-6) 


2.8(-6) 


2.2(-6) 


2.1(-6) 


2.1(-6) 


4.5(-6) 


5.3(-8) 


l.l(-7) 


CO(7-6) 


3.2(-9) 


l.l(-7) 


7.8(-6) 


1.4(-7) 


1.2(-7) 


1.3(-7) 




3.4(-8) 


l.l(-7) 


CO(10-9) 


8.4(-ll) 


1.7(-11) 


2.3(-10) 


9.2(-9) 


4.7(-9) 


2.4(-9) 








CO(16-15) 




1.0(-12) 


1.2(-10) 




4.2(-9) 


7.1(-11) 








13 CO(2-l) 


2.1(-7) 


2.1(-7) 


2.8(-7) 


l.l(-7) 


1.0(-7) 


9.7(-8) 


8.1 (-8) 


1.8(-10) 


1.6(-10) 


[CI] 609 urn 


5.6(-6) 


7.2(-6) 


8.9(-6) 


l.l(-5) 


1.2(-5) 


1.4(-5) 


3.6(-5) 


1.0(-4) 


1.2(-4) 


CO(2-1)/CO(1-0) 


5.9 


5.9 


6.3 


7.5 


7.6 


7.8 


8.4 


25.4 


30.3 


CO(4-3)/CO(l-0) 


17.0 


18.7 


27.2 


26.2 


25.6 


25.8 


36.2 


1.5(2) 


3.6(2) 


CO(7-6)/CO(3-2) 


3.4(-3) 


l.l(-l) 


5.2 


8.6(-2) 


7.9(-2) 


8.0(-2) 




9.8(-l) 


2.2 


CO(10-9)/CO(7-6) 


2.5(-2) 


1.5(-4) 


3.0(-5) 


6.7(-2) 


3.7(-2) 


1.9(-2) 








CO(16-15)/CO(1-0) 




1.2(-5) 


l.l(-3) 




5.0(-2) 


8.7(-4) 








CO(16-15)/CO(10-9) 




6.2(-2) 


0.51 




9.0(-l) 


3.0(-2) 








[CI] 609 Mm/ 13 CO(2-l) 


26.7 


34.1 


31.5 


97.1 


1.2(2) 


1.4(2) 


4.4(2) 


5.6(5) 


7.4(5) 
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Table 3: Line intensities (erg s 1 cm 2 sr 1 ) and ratios at density n = 10 4 cm 



Model 
Radiation field 


Go = 10 s 


Low-CR PDR 

G = 10 4 Go = 10 5 


Go = 10 3 


High-CR PDR 

G = 10 4 Go = 10 s 


F x = 1.6 


XDR 
F x = 16 


F x = 160 


CO(l-O) 


1.6(-7) 


1.9(-7) 


2.5(-7) 


2.1(-7) 


2.4(-7) 


2.7(-7) 


2.6(-7) 


1.4(-10) 


3.4(-ll) 


CO(2-l) 


1.2(-6) 


1.6(-6) 


2.0(-6) 


1.8(-6) 


2.0(-6) 


2.4(-6) 


3.5(-6) 


5.2(-9) 


1.4(-9) 


CO(3-2) 


3.4(-6) 


4.6(-6) 


6.2(-6) 


5.3(-6) 


6.2(-6) 


7.5(-6) 


1.3(-5) 


4.0(-8) 


1.0(-8) 


CO(4-3) 


5.7(-6) 


8.5(-6) 


1.2(-5) 


1.0(-5) 


1.2(-5) 


1.5(-5) 


2.9(-5) 


1.7(-7) 


4.3(-8) 


CO(7-6) 


2.9(-6) 


1.3(-5) 


2.7(-5) 


1.6(-5) 


2.3(-5) 


3.5(-5) 


9.7(-5) 


1.8(-6) 


6.4(-7) 


CO(10-9) 


2.4(-9) 


7.4(-8) 


1.8(-6) 


2.4(-7) 


1.0(-6) 


4.9(-6) 


l.l(-4) 


3.5(-6) 


1.6(-6) 


CO(16-15) 


l.l(-lO) 


6.8(-9) 


2.3(-8) 


2.1(-10) 


3.1(-9) 


2.9(-8) 


l.l(-7) 


2.5(-6) 


1.7(-6) 


HCN(l-O) 


5.2(-9) 


6.8(-9) 


3.8(-9) 


3.9(-9) 


3.8(-9) 


3.8(-9) 


4.9(-10) 


2.2(-13) 


3.1(-14) 


HCN(4-3) 


3.3(-7) 


4.0(-7) 


1.0(-8) 


9.5(-9) 


1.0(-8) 


1.0(-8) 


1.3(-9) 


4.0(-ll) 


3.2(-ll) 


HCO+(1-0) 


1.9(-8) 


2.2(-8) 


2.6(-8) 


l.l(-8) 


1.1 (-8) 


1.1 (-8) 


1.5(-8) 


1.0(-14) 


1.4(-15) 


HCO+(4-3) 


4.3(-7) 


6.6(-7) 


7.9(-7) 


6.2(-8) 


7.0(-8) 


7.2(-8) 


2.0(-7) 


l.O(-ll) 


1.4(-12) 


13 CO(2-l) 


5.8(-7) 


7-4(-7) 


9.3(-7) 


8.0(-7) 


8.8(-7) 


9.6(-7) 


3.0(-7) 


l.O(-lO) 


2.3(-ll) 


[CI] 609 urn 


8.4(-6) 


9.9(-6) 


l.l(-5) 


1.2(-5) 


1.4(-5) 


1.7(-5) 


5.2(-5) 


1.2(-4) 


1.2(-4) 


CO(2-1)/CO(1-0) 


7.8 


8.1 


8.2 


8.4 


8.6 


8.8 


13.3 


36.9 


39.6 


CO(4-3)/CO(l-0) 


35.8 


44.2 


50.4 


47.5 


51.8 


56.9 


1.1(2) 


1.2(3) 


1.3(3) 


CO(7-6)/CO(3-2) 


8.4(-l) 


2.8 


4.3 


3.1 


3.7 


4.7 


7.5 


46.1 


61.9 


CO(10-9)/CO(7-6) 


8.3(-4) 


5.6(-3) 


6.5(-2) 


1.5(-2) 


4.4(-2) 


1.4(-1) 


1.1 


1.9 


2.6 


CO(16-15)/CO(1-0) 


7.0(-4) 


3.5(-2) 


9.2(-2) 


1.0(-3) 


1.3(-2) 


l.l(-l) 


4.1(-1) 


1.8(4) 


5.0(4) 


CO(16-15)/CO(10-9) 


4.6(-2) 


9.2(-2) 


1.2(-2) 


8.8(-4) 


3.1(-3) 


6.0(-3) 


1.0(-3) 


7.0(-l) 


1.0 


HCN(1-0)/CO(1-0) 


3.2(-2) 


3.5(-2) 


3.4(-2) 


1.9(-2) 


1.6(-2) 


1.4(-2) 


1.9(-3) 


1.6(-3) 


9.1(-4) 


HCN(4-3)/CO(4-3) 


5.7(-2) 


4.7(-2) 


3.4(-2) 


9.5(-4) 


8.4(-4) 


6.8(-4) 


4.4(-5) 


2.4(-4) 


7.3(-4) 


HCN(1-0)/HCO+(1-0) 


2.8(-l) 


3.1(-1) 


3.3(-l) 


3.6(-l) 


3.4(-l) 


3.6(-l) 


3.2(-2) 


21.7 


22.6 


HCN(4-3)/HCO+(4-3) 


7.6(-l) 


6.1(-1) 


5.4(-l) 


1.5(-1) 


1.5(-1) 


1.4(-1) 


6.3(-3) 


4.0 


22.6 


[CI] 609 ^m/ 13 CO(2-l) 


14.6 


13.4 


12.1 


14.9 


16.2 


17.6 


1.7(2) 


1.2(6) 


5.1(6) 
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Table 4: Line intensities (erg s 1 cm 2 sr 1 ) and ratios at density n = 10 5 cm 



Model 




Low-CR PDR 






High-CR PDR 






XDR 




Radiation field 


Go = 10 s 


G = 10 4 


Go = 10 5 


Go = 10 3 


Go = 10 4 


Go = 10 5 


F x = 1.6 


F x = 16 


F x = 160 


CO(l-O) 


3.1 (-7) 


4.1 (-7) 


5.5(-7) 


3.3(-7) 


4.3(-7) 


5.6(-7) 


7.6(-7) 


1.4(-6) 


1.2(-6) 


CO(2-l) 


2.7(-6) 


3.5(-6) 


4.7(-6) 


2.8(-6) 


3.7(-6) 


4.9(-6) 


7.1 (-6) 


1.4(-5) 


1.5(-5) 


CO(3-2) 


8.5(-6) 


1.2(-5) 


1.6(-5) 


9.1(-6) 


1.2(-5) 
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Fig. I. — Chemical and thermal structure of PDR and XDR models at density n = 10 5 cm 3 
and G = 10 3 (F x = 1.6 ergs s" 1 cm" 2 ) 



